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In an excitable thin-layer Belousov-Zhabotinsky (BZ) medium a localised perturbation leads to for¬ 
mation of omnidirectional target or spiral waves of excitation. A sub-excitable BZ medium responds 
to asymmetric local perturbation by producing travelling localised excitation wave-fragments, dis¬ 
tant relatives of dissipative solitons. The size and life span of an excitation wave-fragment depend on 
the illumination level of the medium. Under the right conditions the wave-fragments conserve their 
shape and velocity vectors for extended time periods. We interpret the wave-fragments as values of 
Boolean variables. When two or more wave-fragments collide they annihilate or merge into a new 
wave-fragment. States of the logic variables, represented by the wave-fragments, are changed in the 
result of the collision between the wave-fragments. Thus, a logical gate is implemented. Several 
theoretical designs and experimental laboratory implementations of Boolean logic gates have been 
proposed in the past but little has been done cascading the gates into binary arithmetical circuits. 
We propose a unique design of a binary one-bit full adder based on a fusion gate. A fusion gate 
is a two-input three-output logical device which calculates conjunction of the input variables and 
conjunction of one input variable with negation of another input variable. The gate is made of three 
channels: two channels cross each other at an angle, third channel starts at the junction. The chan¬ 
nels contain BZ medium. When two excitation wave-fragments, travelling towards each other along 
input channels, collide at the junction they merge into a single wave-front travelling along the third 
channel. If there is just one wave-front in the input channel, the front continues its propagation 
undisturbed. We make a one-bit full adder by cascading two fusion gates. We show how to cascade 
the adder blocks into a many-bit full adder. We evaluate feasibility of our designs by simulating 
evolution of excitation in the gates and adders using numerical integration of Oregonator equations. 


I. INTRODUCTION 

Information processing potential of a spatially 
extended, thin layer, Belousov-Zhabotinsky (BZ) 
medium mm was firstly demonstrated in late 1980s in 
pioneer works by Kuhnert, Agladze and Krinsky on ex¬ 
perimental laboratory implementation of spatial memory 
devices and image processing prototypes in the excitable 
chemical medium mm- Their findings, albeit with a 
ten years delay, ignited a series of well founded theoret¬ 
ical and experimental works on realisation of computing 
devices in BZ medium. These include logical gates imple¬ 
mented in geometrically constrained BZ medium mm, 
approximation of shortest path by excitation waves EHa, 
memory in BZ micro-emulsion uni, information coding 
with frequency of oscillations m , onboard controllers for 
robots [Hilll, chemical diodes m, BZ neuromorphic ar¬ 
chitectures USHIS], wave-based counters [20], and other 
aspects of information processing in excitable chemical 
systems [2TH24] . 

So far no complete arithmetical circuit has ever been 
implemented in BZ in experimental conditions. Some 
preparatory steps have been done. They include simu¬ 
lation and experimental laboratory realisation of basic 
logical gates [26ll29] and generators of mobile localiza¬ 
tions (they can play a role of constant True) in light- 
sensitive BZ medium [30] and adaptive design of simple 
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logical gates using machine learning techniques [3T] . The 
main reason for such slow experimental progress is insta¬ 
bility of localised excitation wave-fragments: the waves 
are either collapse or expand unless some external peri¬ 
odical stimulation is applied. First one-bit half-adder in 
geometrically-constrained light-sensitive BZ medium is 
proposed in [32], however the implementation required 
dynamical update of the illumination level to prevent 
excitation wave-fragments from collapsing or expanding. 
Models of multi-bit binary adder, decoder and compara¬ 
tor in BZ are proposed in [331136] . These architectures 
employ such crossover structures as T-shaped coincidence 
detectors m and chemical diodes m that heavily rely 
on heterogeneity of geometrically constrained space. 

Indeed by controlling excitability [38] in different loci 
of the medium we can achieve impressive results, as it 
is demonstrated in works related to polymorphic logical 
gates [39], analogs of dendritic trees [19] and, particu¬ 
larly, implementation of four-bit input, two-bit output 
integer square root circuits based on alternating ‘con¬ 
ductivity’ of junctions between channels m, however 
experimental prototypes show a high degree of instability 
and sensitivity to environmental conditions. In present 
paper we try to overcome the problem of excitation 
wave-fragments instability via combining geometrical- 
constraining and collision-based approaches. 

The paper is structured as follows. We define Oreg¬ 
onator model of BZ medium and discuss parameters of 
simulation in Sect. im In Sect, mil we outline differences in 
wave-front behaviour in excitable and sub-excitable BZ 
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medium. The fusion gate is designed and its behaviour 
is simulated and analysed in Sect.|lV| Section [V| presents 
architecture of a one-bit half adder made of a fusion gate. 
Two fusion gates are cascaded into a one-bit full adder 
in Sect. ED In Sect.EID we show how to cascade one-bit 
full adders into a many-bit full adder. Feasibility of the 
approach is discussed in Sect. lYm] 
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II. MODEL 


We use two-variable Oregonator equation m adapted 
to a light-sensitive Belousov-Zhabotinsky (BZ) reaction 
with applied illumination [42] : 


du 

dt 


-(w -u^- ifv + (/>)^) + 
e u-\- q 


dv 

dt 


= u — V 


( 1 ) 


The variables u and v represent local concentrations of an 
activator, or an excitatory component of BZ system, and 
an inhibitor, or a refractory component. Parameter e sets 
up a ratio of time scale of variables u and g is a scaling 
parameter depending on rates of activation/propagation 
and inhibition, / is a stoichiometric coefficient. Constant 
0 is a rate of inhibitor production. In a light-sensitive BZ 
0 represents the rate of inhibitor production proportional 
to intensity of illumination §■ 


We integrate the system using Euler method with five- 
node Laplace operator, time step At = 0.001 and grid 
point spacing Ax = 0.25, e = 0.02 (see Sect. Ill), / = 1.4, 
q = 0.002. The parameter (j) characterises excitability of 
the simulated medium. 


To generate excitation waves of wave-fragments we 
perturb the medium by square solid domains of excita¬ 
tion, 10 X 10 sites in state u = 1.0; if different shape of 
perturbation was used we indicate this in the captions of 
the figures. 

The medium is excitable and exhibits ‘classical’ target 
waves when 0 = 0.05 and the medium is sub-excitable 
with propagating localizations, or wave-fragments, when 
(/) = 0.0766. Time lapse snapshots provided in the paper 
were recorded at every 150 time steps, we display sites 
with u > 0.04. 

The model has been repeatedly verified by us in ex¬ 
perimental laboratory studies of BZ system, and the per¬ 
fect match between the model and the experiments was 
demonstrated in [T7H3Q] . 

We adopt the following symbolic notations. Boolean 
variables x and y take values ‘O’, logical False, and ‘1’, 
logical True; xy is a conjunction (operation and), x^y 
is exclusive disjunction (operation xor), x is a negation 
of variable x (not). Disjunction (OR) of variables x and 
^ is X + 


(b) 


(c) 

FIG. 1. Time lapsed snapshots of wave-fragments propagat¬ 
ing in simulated BZ medium in (a) near lower threshold of 
excitability, 0 = 0.079, (b) sub-excitable mode, (j) — 0.07905, 
(c) non-excitable mode, 0 = 0.0791. The media were per¬ 
turbed by rectangular north-south elongated domains of ex¬ 
citation, 3 X 40 sites in state u — 1.0. Sites of initial segment- 
wise perturbation are shown by arrows. Grid size is 1125x250 
nodes. The pictures are not snapshots of many wave fronts 
generated at the initial stimulation point. These are time 
lapsed snapshots of two waves (one propagating right, an¬ 
other left) recorded every 150^^ step of numerical integration. 


III. EXCITABLE VERSUS SUB-EXCITABLE 
MEDIA 

Given initial asymmetric excitation domain a wave- 
fragment is formed. The fragment’s velocity vector is 
a normal to longest side of the perturbation domain. 
The medium is excited by rectangular domains of per¬ 
turbed sites. The perturbation domains are elongated 
along north-south axis therefore wave-fragments gener¬ 
ated propagate west and east. In excitable medium, 
with low threshold of excitation, every resting site neigh¬ 
bouring excited site becomes excited. Therefore literally 
a point-wise excitation is enough to generate full-scale 
omni-directional wave. In sub-excitable medium a rest¬ 
ing site needs several excited neighbours to get excited. 
Initially asymmetric wave will cause excitation of rest¬ 
ing sites residing nearby middle of the wave-front only. 
Therefore the excitation wave-fragment ‘grows’ from its 
centre and ‘dies’ at its ends. This is why wave-fragments 
in Fig. [^propagate only east and west and shrink or just 
slightly expand along south-north direction. 

Depending on medium’s excitability wave-fragments 
may expand (Fig. keep their shape for a long time 
(Fig. [^) or collapse (Fig.[^). In ideal situation, assum¬ 
ing that wave-fragments keep their shapes indefinitely, we 
can implement a collision-based computing circuit of any 
depth subject to space availability. However in reality, 
particularly in conditions of chemical laboratory exper¬ 
iments, wave-fragments are very unstable. It is almost 
impossible to keep a medium at the precise level of sub- 
excitability (Fig. and almost any wave-fragment will 
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FIG. 2. Time lapsed overlays of the fusion of two excitation 
wave-fragments. One fragment is travelling from north-west 
to south-east, another fragment from south-west to north¬ 
east. The wave-fragments collide and fuse into a new localised 
excitation travelling east. 

expand (Fig. [^) or collapse (Fig. [^). 

It is possible to keep wave-fragment from collapsing 
or exploding by periodically changing excitability of the 
medium [43]. When a wave-fragment expands we in¬ 
crease illumination level thus decreasing the medium’s 
excitability. When the wave-fragment starts to collapse 
the illumination is decreased, the medium’s excitability 
increases and the wave-fragment expands. In designs of 
logical circuits when channels are represented by shad¬ 
owed areas and other parts of medium with illuminated 
areas, it would be difficult to keep periodicity of illumina¬ 
tion and synchronise signals are the same time. There¬ 
fore, we adopted the following approach. We keep the 
medium in excitable mode inside channels and in sub- 
excitable mode at the junctions. Thus we do not worry 
about collapsing of wave-fronts inside channels, yet we 
keep wave-fronts localised when they cross junctions be¬ 
tween channels. 


IV. FUSION GATE 

When two wave-fragments, of equal size, travelling to¬ 
wards each other in a sub-excitable medium collide head- 
on, they annihilate. When the wave-fragments collide 
at an acute angle they fuse into a single wave-fragment 
(Fig.|2|. If one of the wave-fragments was present an¬ 
other fragment would move along its original trajectory. 

This give us an idea of a fusion gate. The gate is 
built of channels containing BZ medium being in a sub- 
excitable mode. The gate has two input channels and 
three outputs channels (Fig. If). When there is a wave- 
fragment present in the input channel x (y) we assume 
that input variable x (y) takes logical Boolean value 
True, x = 1 {y = 1); otherwise Boolean value False, 
X = 0 {y = 0). 

If the gate is in excitable mode and initiated with in¬ 
put X = 1 and ^ = 0 (Fig. If and Fig. |^b), then wave- 
fragment propagating along input channel x expands at 
the junction (Fig. |^d). Excitation wave-front then en¬ 
ters all other channels (Fig. If and Fig. Iff). Such gate 
implements only multiplication of a signal and it is not 
useful for a realisation of logical functions. If the gate 
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FIG. 3. Fusion gate, (a) Scheme: inputs are x and y, out¬ 
puts are xy, xy, xy. (a-d) Time lapsed overlays of excitation 
waves, (b) Excitable mode, x = 1, y = 0. (c) Sub-excitable 
mode, X = 1, y = 0. (d ) Excitable mode, x = 1, y = 1. 
(e) Sub-excitable mode, x = 1, y = 1. 


is in a sub-excitable mode and initiated with x = 1 and 
y = 0 (Fig. |3f and Eig. Ifh) then the wave-fragment 
does not expand at the junction (Eig. |^j) but contin¬ 
ues through the junction along its original trajectory 

Fig. If 1). 

If the gate was in excitable mode and initiated with 
input X = 1 and y = 1 (Eig. If) the wave-fragments 
collide, fuse and expand. The excitation wave-fronts en¬ 
ter all three output channels. In sub excitable mode 
the wave-fragments fuse (Eig. |^) into a single localised 
wave-fragment which enters the central output channel 

(Fig.lJ). 

Thus, assuming that presence of a wave-fragment in 
an output channel symbolises logical variable in state 
True at the output channel, we conclude that the fusion 
gate, being in sub-excitable mode, implements two-input- 
three-output Boolean logic gate {x, y) l^y, xy, xy) 
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(j)t = 10.135 (k)t = 11.31 (l)t = 13.605 


FIG. 4. Snapshots of fusion gate for inputs x = 1 and y = 0. (a-f) Excitable mode, (g-1) Sub-excitable mode. Time is shown 
in captions. Excitable heads of the wave-fragments are shown in red. Refractory tail in blue. 


(Fig#). The fusion gate is a key component of a one-bit V. HALF-ADDER 

half-adder presented in next section. 

A one-bit half-adder is a device with two inputs x and 
y and two outputs xy (Carry out) and x ^ y (Sum). 
The half-adder based on the fusion gate is shown in 
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(d)t = 0.475 (e)t = 0.475 


FIG. 5. Snapshots of fusion gate for inputs x = 1 and y — 1. The system is in sub-excitable mode. Time is shown in captions. 
Excitable heads of the wave-fragments are shown in red. Refractory tail in blue. 



x©y xy 


FIG. 6. A scheme of a one-bit half-adder: input channels are 
a and 5, output channels are g and h, internal channels d, e, 
/, and junctions are c, i and j. Input variables x and y are 
fed into channels a and 5, results x ®y and xy are read from 
channels g and h. 


Fig-i It consists of two input channels a and b and 
two output channels g and h. Presence/absence of a 
wave-fragment in an input/output channel symbolises 
logical True/False state of the input variable assigned 
to the channel. Synchronisation of signal wave-fragments 
is achieved geometrically: \a\ + |e| + \g\ = |6| + \d\ + \g\ = 



(c) (d) 


FIG. 7. Time lapsed overlays of excitation waves propagat¬ 
ing through junctions i (ab) and j (cd). (ab) Wave-fragment 
approaches junction i via channel / (Fig. 1^. (cd) Wave- 
fragment approaches junction j via channel d (Fig. [^. 
(ac) The system is in excitable mode, (bd) The system is 
in sub-excitable mode. 


\a\ + I/I + |/i|, where | • | is a length of a channel (Fig.[^. 
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The half-adder has three junctions c, i and j (Fig.J^. 
The junction c is key part of the fusion gate (Fig. [^) 
and propagation of wave-fragments through this junction 
has been discussed in Sect. Propagation of waves via 
junctions i and j is illustrated in Fig. for excitable 
(Fig#c) and sub-excitable (Fig. [^d) modes. 

If the system was in excitable mode then wave- 
fragment entering junction i via channel / would expand 
and propagate into channels e (above and below the junc¬ 
tion) and h (Fig. 1^). Similarly, the wave-fragment en¬ 
tering junction j from channel d expands and propagates 
into channels j, e and h (Fig. Signals are multiplied 
and thus, again, the device being in excitable mode is not 
useful for implementing logical functions. When the sys¬ 
tem is in sub-excitable mode a wave-fragment entering 
junction i via channel / stays localised and propagates 
only into channel h (Fig. Similarly, a wave-fragment 
entering junction j from channel d propagates only into 
channel g (Fig. W)- There is no such combination of in¬ 
puts which might lead to collision of two wave-fragments 
in junction i, therefore we do not consider such scenario. 

The half-adder in action is shown in Fig. When in¬ 
puts are x = 1 and ^ = 0 a wave-fragment is initiated 
in channel a. The wave stays localised (without expan¬ 
sion) when propagating via junction c and enters chan¬ 
nel e. The wave travels along the channel e and through 
junction i and j and enters output channel g (Fig. [^). 
Similarly, when inputs are x = 0 and ^ = 1 a wave is ini¬ 
tiated in channel 6, propagates though junction c, along 
channel d, via junction j and enters output channel g 
(Fig. [^). If both inputs takes value True excitation 
wave-fragments are initiated in channels a and b simul¬ 
taneously. The wave-fragments collide at junction c and 
fuse into a single wave-fragment. This newly born wave- 
fragment propagates into channel /. The wave travels 
along channel /, passes though junction i without ex¬ 
pansion and enters the output channel h (Fig. |^). 

A wave-fragment appears in output channel g only if 
either a wave-fragment was initiated either in input chan¬ 
nel a or input channel b. Therefore, the channel g repre¬ 
sents logical operation x ^ y. A wave-fragment appears 
in output channel h only if wave-fragments were initiated 
in input channels x and y. Therefore, the channel h rep¬ 
resents logical operation xy. Thus, we illustrated how 
the half-adder shown in Fig. is derived. 


VI. ONE-BIT FULL ADDER 

A one-bit full adder is a device with three inputs x, y 
and z (Carry in) and two outputs xyAz{x^y) (carry out) 
and x^y^z (Sum). The full adder is made of two half¬ 
adders (Fig.[^ by cascading Carry in into input channel 
of second half-adder. Sum of first half-adder into second 
input of second half-adder, and Carry out of first half¬ 
adder into OR junction with Carry out output channel 
of second half-adder. The one-bit full adder has three 
inputs channels a, b and c (they represent carry in and 


TABLE I. Schematic dynamic of one-bit full adder for all 
possible combinations of inputs. The adder, as shown in Fig.|^ 
is rotated to the left. 


X 

y 

Cin 

Sum 

Cout 

Configuration 

0 

0 

0 

0 
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1 

0 

0 

1 
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0 
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0 
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0 

1 

1 
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1 

0 

1 

0 

1 








0 

1 

1 

0 

1 








1 

1 

1 

1 

1 



added bits x and y) and two output channels t and u 
(they represent sum x ^y ^ z and carry out xy A z{x ^ 
y)). Synchronisation of signals is achieved geometrically: 
\a\ -h |e| + \p\ A \t\ = |6| -1- \h\ A \k\ + \n\ + |t| = |c| + |/| + 
\k\ + |n| + |t| = |6| -h 1^1 + |/| + \u\. Six junctions d, i, 
j, m, r and s are types of junctions already discussed in 
Sect. El 


The one-bit full adder in action is shown in Fig. 10 


When Carry in has bit up, 2 : = 1, and two inputs bits are 
down, X = 0 and y = z the wave-fragment is originated 
only in input channel a (Fig. pT^). The wave travels 
along channel e, through junction m, enters channel p, 
propagates though junction r and finishes its journey in 
output channel t. For input tuple z = 0^ x = 1 and p = 0, 
an excitation wave-front is initiated in input channel 6, 
propagates through junction d into channel h, through 
junctions j and i into channel k, through junction n into 
channel n, and finally enters the output channel t. When 
inputs are z = 0,, x = 0 and p = 1 a wave-front initiated 
in input channel c propagates through d, along /, through 
i, along k, through m, along n and ends up in output 
channel t. For combination of inputs z = 0^ x = 1 and 
y = I wave-fragments are initiated in input channels b 
and c, they merge into a single wave-fragment at junction 
d. This newly born wave-fragment propagates along p. 
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(a) (b) (c) 


FIG. 8. Time lapsed overlays of excitation waves propagation in the one-bit half-adder (Fig. for inputs (a) x = 1, y = 0, 
(b) X = 1, y = 0, (c ) X = 1, y = 1. Sites of initial segment-wise perturbation are visible as discs. Grid size is 500x790 nodes. 
The pictures are not snapshot of many wave fronts generated at the initial stimulation point. These are time lapsed snapshots 
of a single wave (ab) or two waves merging into a single wave (c) recorded every 150^^ step of numerical integration. 


through along /, through s and into output channel u. 
Dynamics of excitation waves for inputs z = z = 1 and 
^ = 0 is shown in Fig. pTj p. The wave-fragments merge 
into a single wave-fragment at junction m and the newly 
born fragment propagates into output channel u. When 
all three input has state bit up, z = x = y = the 
wave-fragment representing z = 1 does not interact with 
the wave-fragment produced by merged wave-fragments 
representing x = 1 and y = I (Fig. [T^). Thus we have 
wave-fragments appearing on both output channels t and 
u. 

Trajectories of wave-fragments for all possible combi¬ 
nations of inputs are shown in Tab. [Ij The output chan¬ 
nel t symbolises Sum x 0 ^ 0 z and channel u Carry out 
xy 0 z{x 0 y). 


VII. CASCADING ONE-BIT FULL ADDER IN 
MANY-BITS FULL ADDER 

The full adder blocks (Fig.|^ can be straightforwardly 
cascaded into a many-bit adder. An example of a 4- 
bit adder is shown in Fig. As we can see, an n-bit 
full adder can be made of 2n one-bit half-adders. Sig¬ 
nals, represented by wave-fragments, are fully synchro¬ 
nised. Therefore, we can feed input data into the many- 
bit adder with period exceeding two wave-width, so ex¬ 


cited heads of current signals do not interfere with re¬ 
fractory tails of previous signals. Assuming that a width 
of excitation wave-front in BZ system is circa 0.5-1 mm 
m ns, and the speed of propagation is about 1 mm 
per minute, we believe that many-bits full adder’s clock 
frequency could be 0.01 Hz. With regards to a spatial 
complexity, to stay on the safe side we assume that one- 
bit full adder fits into 10 x 3 mm box, thus n-bit adder 
would take 10 • n x 3 mm space. That is e.g. eight-bit 
adder implemented with fusion gates in sub-excitable BZ 
medium would be 8 cm by 0.3 cm in size. 


VIII. DISCUSSION 

We proposed implementation of binary adders in 
architectures of homogeneous channels in Belousov- 
Zhabotinsky (BZ) system. A fusion gate is a key el¬ 
ement of the architecture. The function of the fusion 
gate is to allow lonely excitation wave-fragments to pass 
along their original trajectories undisturbed but direct¬ 
ing a new wave-fragment born after collision between 
two wave-fragments along the new trajectory. Localised 
wave-fragments exist only in sub-excitable BZ system. 
Unfortunately, they are unstable: after some period of 
propagation the wave-fragments start expanding or col¬ 
lapsing. Thus to keep signals propagation without de- 


z X _y 



x@y@z xy+z(x@y) 


FIG. 9. A scheme of a one-bit full-adder: input channels are 
a, b and c, output channels are t and internal channels e, 
/, h, /c, /, n, o, p, and junctions are d, z, j, m, r, s. Input 
variables a: and p are fed in channels h and c, Carry In is fed 
in channel a, results x®y®z and xy-{-z[x®y) are read from 
channels t and u. 

cay, we kept BZ medium in the excitable mode inside 
the channel. To prevent signals from exploding at the 
junction we kept the BZ medium in sub-excitable mode 
at the junctions between the channels. If required, the 
whole architecture can be kept in sub-excitable mode by 
enforcing periodic changes of excitability. 

In experimental laboratory conditions, the architecture 
can be implemented in two ways. We can make physical 
templates of the arithmetical circuits, where channels are 
depressions filled with BZ solution. Such approach would 
be robust, however the whole architecture would be non- 
reconfigurable. And we would still need to apply illumi¬ 
nation to control excitability. Using only light would be 
the most efficient and elegant approach. If BZ medium is 
light-sensitive then generation and propagation of waves 
is controlled by illlumination. In the Oregonator model 
0 this control is expressed via constant (j). This is a 
rate of inhibitor production proportional to intensity of 
illumination. Thus the architecture of a computing cir¬ 
cuit can be projected onto the medium, in such a manner 
that channels are black and background is white. Then 


excitation waves will be confined to the channels because 
rest of the medium would be non-excitable. Moreover, by 
dynamically varying illumination at the junctions we can 
implement polymorphic logical gates by controlling out¬ 
comes of inter-fragment collisions using the illumination 
level. For example, in [39] we demonstrated that how 
to switch the BZ collision gate between NOR to XNOR 
modes. 

Let us compare our design with previous implemen¬ 
tations (Fig. Our previous design of one-bit half¬ 
adder implemented in BZ system in computer model [3^ 
and laboratory experiments [46| has four inputs (data 
signals are duplicated) and four junctions where wave- 
fragments interact (Fig . [T^ ). Half-adder implemented 
in computer models in |35j has two inputs and two out¬ 
puts and seven junctions, or interactions loci (Fig. [T^). 
The architecture [35| consists of six disconnected domains 
of ‘conductivity’ and has seven junctions, which control 
propagation of excitation wave-fronts. The half-adder 
designed in our present paper (Fig. [T^ ) has two input 
and two outputs, only three junctions between channels, 
and wave-fragment can interact only in two junctions. 

Obviously, one can find advantages in every implemen¬ 
tation. For example, the first ever BZ half-adder de¬ 
sign [32j [46] produces conjunction in addition to carry 
out and sum; and, the design [35| does not require a 
control of excitability at the junctions. Yet our present 
architecture seems to be more optimal because it is made 
of a single domain of ‘conductivity’ where excitation 
wave-fragments can interact only in two junctions. Also, 
speaking in terms of laboratory experiments each coin¬ 
cidence detector [37] and/or chemical diode [15], used 
in [35], add to overall unreliability of the system be¬ 
cause things can go wrong when excitation waves try to 
pass loci of discontinuity, because even a width of chan¬ 
nel [47| [48] strongly affects dynamics of wave-fronts. 

There are three obvious directions of further studies: 
speed up, application domain, and extension to many¬ 
valued circuits. BZ circuits implemented in experimental 
laboratory conditions are slow, each cycle of computation 
takes quarter of an hour. To speed up the computation, 
we can implement the BZ adders in analog devices [53] . 
molecular arrays m and at the nano-scale [54]. With 
regards to the application domain, we envisage that our 
designs can be successfully integrated in embedded chem¬ 
ical sensing devices m, m- And yet another direc¬ 
tion of future research could be in expanding our designs 
to many-valued logic circuits. In principle, we can im¬ 
plement three-valued gates with wave-fragments in BZ 
system [58], based on facilitation or inhibition between 
two subsequent excitation fronts, more work is required 
though to achieve cascading of these many-valued logic 
gates. 
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FIG. 10. Time lapsed overlays of excitation waves propagation in the one-bit full adder (Fig. for inputs (a) x = 0, y = 0, 
z = 1, (b) X = 1, y = 0, z = 1, (c ) X = 1, y = 1., z = 1. 
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